A nonextensive approach to the dynamics of financial observables 
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We present results about financial market observables, specifically returns and traded volumes. 
They are obtained within the current nonextensive statistical mechanical framework based on the 



entropy S q — k- 



£ Pi w 
1=1 {q G R) (Si = Sbg = ~kJ2 Pl ki Pi 



1-9 



More precisely, we present stochas- 



tic dynamical mechanisms which mimic probability density functions empirically observed. These 
mechanisms provide possible interpretations for the emergence of the entropic indices q in the time 
evolution of the corresponding observables. In addition to this, through multi-fractal analysis of 
return time series, we verify that the dual relation q s tat + Isens = 2 is numerically satisfied, q B tat 
and q se ns being associated to the probability density function and to the sensitivity to initial con- 
ditions respectively. This type of simple relation, whose understanding remains ellusive, has been 
empirically verified in various other systems. 



I. INTRODUCTION 

In recent years statistical mechanics has enlarged its 
original assignment: the application of statistics to large 
systems whose states are governed by some Hamilto- 
nian functional |lj . Its capability for relating microscopic 
states of individual constituents of a system to its macro- 
scopic properties are nowadays used ubiquitously 
Certainly, the most important of these connections still is 
the determination of thermodynamic properties through 
the correspondence between the entropy concept, origi- 
nally introduced by Rudolf Julius Emmanuel Clausius in 
1865 Pi, and the number of allowed microscopic states, 
introduced by Ludwig Boltzmann around 1877 when he 
was studying the approach to equilibrium of an ideal 
gas 0. This connection can be expressed as 



S = k lnW, 



(1) 



where k is a positive constant, and W is the number 
of microstates compatible with the macroscopic state of 
an isolated system. This equation, known as Boltzmann 
principle, is one of the cornerstones of standard statistical 
mechanics. 

When the system is not isolated, but instead in con- 
tact to some large reservoir, it is possible to extend Eq. 
i[TJl. under some assumptions, and obtain the Boltzmann- 
Gibbs entropy 



'BG 



w 

i=i 



(2) 
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where nt is the probability of the microscopic configura- 
tion i Jy. The Boltzmann principle should be derivable 
from microscopic dynamics, since it refers to microscopic 
states, but the implementation of such calculation has 
not been yet achieved. So, Boltzmann-Gibbs (BG) sta- 
tistical mechanics is still based on hypothesis such as the 
molecular chaos |4| and ergodicity [5j . In spite of the lack 
of an actual fundamental derivation, BG statistics has 
been undoubtedly successful in the treatment of systems 
in which short spatio/temporal interactions dominate. 
For such cases, ergodicity and (quasi-) independence are 
favoured and Khinchin's approach to Sbg is valid 
Therefore, it is entirely feasible that other physical en- 
tropies, in addition to the BG one, can be defined in 
order to properly treat anomalous systems, for which the 
simplifying hypothesis of ergodicity and/or independence 
are not fulfilled. Examples are: metastable states in 
long-range interacting Hamiltonian dynamics, metaequi- 
librium states in small systems (i.e., systems whose num- 
ber of particles is much smaller than Avogrado's num- 
ber), glassy systems, some types of dissipative dynamics, 
and other systems that in some way violate ergodicity. 
This includes systems with non-Markovian memory (i.e., 
long-range memory), like it seems to be the case of finan- 
cial ones. Generically speaking, systems that may have 
a multi-fractal, scale-free or hierarchical structure in the 
occupancy of their phase space. 

Inspired by this kind of systems it was proposed in 
1988 the entropy @ 



w 

i=l 
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(q e K) , 



(3) 



which generalises Sbg (lim g ^i S q — Sbg), as the basis of 
a possible generalisation of BG statistical mechanics 0,0- 
The value of the entropic index q for a specific system 
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is to be determined a priori from microscopic dynamics. 
Just like Sbgi S q is nonnegative, concave, experimentally 
robust (or Lesche-stable (Vg>0), and leads to a 

finite entropy production per unit time HEl. Mor 

eover, 

it has been recently shown [13 that it is also extensive, 



i.e. 



N 



S q (A 1 + A 2 + ... + A N )~J2 S i ( A >) . 



(4) 



(x - i± q y 



$[p{x)} q dx 



dx = ((x - fj, q f^ = a 2 q , (8) 



corresponding to the generalised mean and variance of x, 
respectively [tj- 

From the variational problem using (j3J) under the 
above constraints, we obtain 



for special kinds of correlated systems, more precisely 
when the phase-space is occupied in a scale-invariant 
form. By being extensive, for an appropriate value of 
q, S q complies with Clausius' concept on macroscopic 
entropy, and with thermodynamics. 

Since its proposal, entropy has been the source of 
several results in both fundamental and applied physics, 
as well as in other scientific areas such as biology, chem- 
istry, economics, geophysics and medicine 01 . Herein, 
we both review and present some new results concern- 
ing applications to the dynamics of financial market ob- 
servables, namely the price fluctuations and traded vol- 
umes. Specifically, we will introduce stochastic dynam- 
ical mechanisms which are able to reproduce some fea- 
tures of quantities such as the probability density func- 
tions (PDFs) and the Kramer-Moyal moments. More- 
over, we will present some results concerning the return 
multi-fractal structure, and its relation to sensitivity to 
initial conditions. 

Our dynamical proposals will be faced to empirical 
analysis of 1 minute returns and traded volumes of the 
30 companies that were used to compose the Dow Jones 
Industrial Average (DJ30) between the 1 st July and the 
31 st December 2004. In order to eliminate specious be- 
haviours we have removed the well-known intra-day pat- 
tern following a standard procedure ||. After that, the 
return values were subtracted from its average value and 
expressed in standard deviation units, whereas the traded 
volumes are expressed in mean traded volume units. 



II. VARIATIONAL PRINCIPLE USING THE 
ENTROPY S q 

Before dealing with specific financial problems, let us 
analyse the probability density function which emerges 
when the variational principle is applied to S q . 

Let us consider its continuous version, i.e., 



S a — h 



1-J[p(x)] q dx 



(5) 



The natural constraints in the maximisation of © are 
p (x) dx = 1, (6) 



(7) 



corresponding to normalisation, and 

[p(x)r 



l\p{x)] q dx 



dx = (x) = n q , 



p(x) = A q 1 + (q - 1) B q (x - fi q y 
where, 



(<?<3) : 



An 
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«<1 

q>l 



and 



B q =[(3-q) 5*\ 



(9) 



(10) 



(11) 



Standard and generalised variances, a 2 and a q respec- 
tively, are related by 



3q 



(12) 



3-q 

Defining the q- exponential function as 

e' q = [l + (l- q )x]^ (e? = e»)> (13) 
(e| = if 1 + (1 — q)x < 0) we can rewrite PDF I© as 

P (x)=A q e- B ^-^ 2 , (14) 

hereafter referred to as q- Gaussian. 

For q — tt 211 , the g-Gaussian form recovers the Stu- 

^ l+m ' * 

dent's ^-distribution with m degrees of freedom (m = 
1,2,3,...) with finite moment up to order m th . So, 
for q > 1, PDF (|14fl presents an asymptotic power- 
law behaviour. On the other hand, if q — with 
n = 3,4,5,..., p(x) recovers the r-distribution with n 
degrees of freedom. Consistently, for q < 1, p(x) has 
a compact support which is defined by the condition 

\x - p q \ < 



l-q °<T 



III. APPLICATION TO MACROSCOPIC 
OBSERVABLES 

A. Model for price changes 

The Gaussian distribution, recovered in the limit 9 — > 1 
of expression i|14fl . can be derived from various stand- 
points. Besides the variational principle, it has been de- 
rived, through dynamical arguments, by L. Bachelier in 
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his 1900 work on price changes in Paris stock market |14j . 
and also by A. Einstein in his 1905 article on Brownian 
motion [15J. In particular, starting from a Langevin dy- 
namics, we are able to write the corresponding Fokker- 
Planck equation and, from it, to obtain as solution the 
Gaussian distribution. Analogously, it is also possible, 
from certain classes of stochastic differential equations 
and their associated Fokker-Planck equations, to obtain 
the distribution given by Eq. lfH|) . 

In this section, we will discuss a dynamical mechanism 
for returns, r, which is based on a Langevin-like equation 
that leads to a PDF (^-Gaussian) with asymptotic power- 
law behaviour 0, 01 • This equation is expressed as 



dr = —krdt ■ 



^0 \p(r,t)f- q) dW t (g>l) : 



(15) 



(in Ito convention) where Wt is a regular Wiener process 
and p{r, t) is the intantaneous return PDF. In a return 
context the deterministic term of cq. H15|) intends to rep- 
resent internal mechanisms which tend to keep the mar- 
ket in some average return or, in a analogous interpreta- 
tion, can be related to the eternal competition between 
speculative price and the actual worth of a company. In 
our case, we use the simplest approach and write it as 
a restoring force, with a constant fc, similar to the vis- 
cous force in the regular Langevin equation. In regard to 
the stochastic term, it aims to reproduce the microscopic 
response of the system to the return: is the volatility 
constant (intimately associated to the variance of p(r, t)) 
and q, the nonextensive index, reflects the magnitude of 
that response. Since the largest unstabilities in the mar- 
ket are introduced by the most unexpected return values, 
it is plausible that the stochastic term in Eq. H15|) can 
have such inverse dependence on the PDF p(r,t). Fur- 
thermore, Eq. (|15|l presents a dynamical multiplicative 
noise structure given by, 

r(t)= [ e-^'-'V^M)]^ dWf, (16) 

J — oo 

where we have assumed r (— oo) = 0. 

The associated Fokker-Planck equation to Eq. (|15fl is 
given by 



dp(r, t) d 



^rp{r,t)] + \^- 2 {e[p(r,t)f^}, 

(17) 



at 



and the long-term probability density function is|lq. \li 



p 0) = -| 



i-a-g) 



kr 2 



(2 - q) Z1- 



(18) 



One of the most interesting features of eq. fT5|l is its 
aptitude to reproduce the celebrated U-shape of the 2 nd 
(i.e., n = 2) Kramers-Moyal moment 



It is this fact which allowed the establishment of analogies 
(currently used in financial mimicry) between financial 
markets dynamics and fluid turbulence |20j |. 

It is noteworthy that eq. (|l(j[) is statistically equivalent 
to [H |2]| 



dr 



-krdt + A q (t) dWt + {q-l) B q (r, t) dW[ , (20) 



i.e., a stochastic differential equation with independent 
additive and multiplicative noises. If eq. I|15|l allows an 
immediate heuristic relation between q and the response 
of the system to its own dynamics, eq. (|20|l permits 
a straighforward dynamical relation between q and the 
magnitude of multiplicative noise in such a way that, for 
q = 1, the Langevin equation is recovered as well as the 
Gaussian distribution. 

In Fig. Q we present the typical PDF for the 1 minute 
returns of a company constituent of the Dow Jones Indus- 
trial Average 30 (upper panel) presenting q = 1.31 ±0.02, 
a time series generated by eq. (|15l) (middle panel), and 
the U-shaped 2 nd Kramers-Moyal moment for our data 
(lower panel). As it can be seen the accordance using the 
simplest approach is already quite nice. Upgrades of this 
model can be obtained by taking into account the risk- 
aversion effects, which induce asymmetry on the PDF, 
and correlations on the volatility in a way which differs 
from others previously proposed. The formulation pre- 
sented herein has also the advantage of being aplicable 
to systems which are not in a stationary state since the 
time-dependent solutions of the Fokker-Planck equation 
are of the g-Gaussian type as well. 



B. Model for traded volumes 

Changes in the price of a certain equity are naturally 
dependent on transactions of that equity and thus on its 
traded volume, v. Previous studies proved the asymp- 
totic power-law behaviour of traded volume PDF [22|, 
later extended for all values of v • In this case it was 
shown that the traded volume PDF is very well described 
by the following ansatz distribution 



P(v) 



1 

Z 



exp 



(21) 



M n (r, t, t) = / (/ - r) n P (/, t + r|r, t) dr' 



[P (»•,*)] 
(19) 



where v represents the traded volume expressed in its 
mean value unit (V), i.e., v = V/(V), p and ip are pa- 
rameters, and Z — J °° exp ? ^— ^\ dv. 

The probability density function l|21|) was recently ob- 
tained from a mesoscopic dynamical scenario [24[ based 
in the following multiplicative noise stochastic differential 
equation 

dv = -7(1; - -)dt + J 2— vdW t , (22) 
a \ a 

(.1-1) 

where W t is a regular Wiener process following a nor- 
mal distribution, and v > 0. The right-hand side terms 
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FIG. 1: Upper panel: Probability density function vs. r. 
Symbols correspond to an average over the 30 equities used 
to built DJ30 and the line represents the PDF obtained from 
a time series generated by eq. 1161 which is presented on 
middle panel. Lower panel: 2 nd Kramers-Moyal moment 
Mi » t0 [p{r)} (1 ~ q) = r^r q [(5 - 3 q) a 2 + (q - 1) r 2 ] from 
which k parameter is obtained and where the stationary hy- 
pothesis is assumed (to = —oo <C — fc _1 <C 0). Parameter 
values: r = 1 min, k = 2.40 ± 0.04, a = 0.930 ± 0.08 and 
q = 1.31 ±0.02. 



of eq. I|22ll represent inherent mechanisms of the sys- 
tem in order to keep v close to some "normal" value, 
Lu/a, and to mimic microscopic effects on the evolution 
of v, like a multiplicative noise commonly used in inter- 
mittent processes. This dynamics, and the correspond- 
ing Fokker-Planck equation [l^, lead to the following 




1000 



FIG. 2: Symbols represent the average correlation function 
for the 30 time series analysed and the line represents a double 
exponential fit with characteristic times of T\ — 13 and T2 = 
332 yielding a ratio about 25 between the two time scales 



Eq. J2H (R 2 

in minutes). 



= 0.991, x = 9 x 10 , and time is expressed 



inverted Gamma stationary distribution: 



/(«) = 



1 



loT [a + 1] 



-a-2 



exp 



(23) 



Consider now, that instead of being a constant, u> is a 
time dependent quantity which evolves on a time scale 
T larger than the time scale of order 7 _1 required by 
eq. H22f> to reach stationarity [2M [2(| . This time depen- 
dence is, in the present model, associated to changes in 
the volume of activity (number of traders that performed 
transactions) and empirically justified through the anal- 
ysis of the self-correlation function for returns. In Fig. |3 
we have verified that the correlation function is very well 
described by 



C [v (t) ,v(t + r)} = Ci e- r / Tl + C 2 e~ T / T2 



(24) 



with T2 = 332 3> T\ = 13. In other words, there is 
first a fast decay of C [v (t) , v (t + r)], related to local 
equilibrium, and then a much slower decay for larger r. 
This constitutes a necessary condition for the application 
of a superstatistical model [2f| . 

If we assume that u> follows a Gamma PDF, i.e., 



P(w) 



1 



XT[S] 



/U}\ S -~ L 




(x) exp 





(25) 



then, the long-term distribution of v will be given by 
p (v) — J / (v) P (oj) du. This results in 
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(i) 


exp ? 








V 



(26) 



where A = (q — 1), 6 = -^j — a — 1. Bearing in mind 
that, for q > 1, 



x 1- 1 e 



(27) 
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FIG. 3: Upper panel: Excerpt of the time series generated 
by our dynamical mechanism (simulation) to replicate 1 min 
traded volume of Citigroup stocks at NYSE (data). Lower 
panel: 1 min traded volume of Citigroup stocks probability 
density function vs. traded volume. Symbols are for data, and 
solid line for the replica. Parameter values: 9 = 0.2 12 ±0.003, 
p = 1.35 ± 0.02, and q = 1.15 ± 0.02 ( X 2 = 3.6 x 10" 4 , 
R 2 = 0.994). 



we can redefine our parameters and obtain the q-Gamma 

PDF (23). 

In Fig. |3| we present a comparation between the traded 
volume of Citigroup (2004 world's number one company 
|27j ) stocks, as well as a replica of that time series ob- 
tained using this dynamical proposal. As it can be easily 
verified, the agreement is remarkable. 

IV. THE NONEXTENSIVE g-TRIPLET AND 
FINANCIAL OBSERVABLES 



Analogously, it was recently conjectured that, for 
systems which can be studied within nonextensive statis- 
tical mechanics, the energy probability density function 
(associated to stationarity or (meta) equilibrium), the 
sensitivity to the initial conditions, and the relaxation 
would be described by three entropic indices q s tat, qsens, 
and q re i, referred to as the q-triplet. The first physi- 
cal corroboration of such scenario has been made from 
the analysis of two sets of daily averages of the mag- 
netic field strength observed by Voyager 1 in the solar 
wind [2!| . Others systems are currently on study (e.g., 
|3fj|). Of course, if the system is non Hamiltonian, it has 
no energy distribution, hence q s tat cannot defined in this 
manner. We may however estimate it through a station- 
ary state generalised Gaussian (which would generalise 
the Maxwellian distribution of velocities for a BG sys- 
tem in thermal equilibrium). In contrast, the other two 
indices, q sen s and q re i, remain defined in the usual way. 

Let us focus now on the multi-fractal structure of re- 
turn time series. It has been first conjectured, and later 
proved, for a variety of nonextensive one-dimensional sys- 
tems, that the following relation holds |3l| : 



1 qsens hxxlin ^max 

where /i m in and h max are respectively the minimal and 
maximal /i-values of the associated multifractal spec- 
trum f(h). In fig. 01 we depict the multifractal spec- 
trum of 1 minute traded volumes, obtained by the appli- 
cation of the MF-DFA5 method H2; h and / (h) have 
been obtained from averages of the empirical data of 30 
companies. Through this analysis, we have determined 
h min = 0.28 ± 0.04 and h max = 0.83 ± 0.04. The use 
of Eq. yields q sens = 0.58 ± 0.10. Considering that 
the q value obtained for the return probability density 
function was q s t a t = 1-31 ± 0.02, we verify that the dual 
relation 

qstat T" qsens — 2 (29) 

is approximately satisfied within the error intervals. Tak- 
ing into account the well-known fast decay of return self- 
correlations, we see that the price changes for a typical 
DJ30 stock may be essentially described by the q-triplet 

{qsen, qstat, qrel} = {0.58 ± 0.10, 1.31 ± 0.02, 1}. 



Systems characterised by Boltzmann-Gibbs statistical 
mechanics present the following characteristics: (i) Their 
PDF for energies is proportional to an exponential func- 
tion in the presence of a thermostat; (ii) They have strong 
sensitivity to the initial conditions, i.e., this quantity in- 
creases exponentially with time (currently referred to as 
strong chaos), being characterised by a positive maxi- 
mum Lyapunov exponent; (iii) They typically present, 
for basic macroscopic quantities, an exponential decay 
with some relaxation time. In other words, these three 
behaviours exhibit exponential functions (i.e., q = 1). 



V. FINAL REMARKS 

In this article we have presented a nonextensive sta- 
tistical mechanics approach to the dynamics of finan- 
cial markets observables, specifically the return and the 
traded volume. With this approach we have been able 
to present mesoscopic dynamical interpretations for the 
emergence of the entropic index q frequently obtained by 
a numerical adjustment for data PDF of eqs. i|18|l and 
1)21(1. For the case of returns, q is related to the reaction 
degree of the agents on the market to fluctuations of the 
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FIG. 4: Multi-fractal spectrum / (h) vs. h for 1 minute return 
averaged over the 30 equities with hmin = 0.28 ± 0.04 and 
U = 0.83 ± 0.04. 



observable, while for the case of traded volume it is asso- 
ciated to fluctuations on the (local) average traded vol- 
ume. Along with these dynamical scenarios, and based 
on the multi-fractal nature of returns, we have verified 
that this quantity appears to approximatively satisfy the 
dual relation, q s tat + Qsens — 2, previously conjectured 
within the emergence of the q-triplet which characterises 



the stationary state, the sensitivity to initial conditions, 
and the relaxation for nonextensive systems. The com- 
plete understanding of these connections remains ellu- 
sive. For instance, concerning relaxation and the q-triplet 
conjecture, a new question arise for price changes. It is 
well-known that the self-correlation for returns is of ex- 
ponential kind, in contrast with the long-lastin g co rrela- 
tions for the volatility (or returns magnitude) 33]. The 
latter is also considered a stylised fact and it is compati- 
ble with a g-exponential form. In this way, if the efficient 
market hypothesis is considered the key element in finan- 
cial markets, then it makes sense to assume q re i — 1. But, 
if arbitrage on markets is considered as the fundamental 
feature instead, then the essential relaxation to be taken 
into account might be the one related to the volatility, for 
which q re i > 1. Progress is clearly still needed, at both 
the fundamental and applied levels, in order to achieve a 
deep understanding of this complex system. 
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